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1 Abstract 

Cold neutrons are converted to ultra-cold neutrons (UCN) by the excitation of a single phonon or multiphonons in superfluid 
helium. The dynamic scattering function S{q, huj) of the superfluid helium strongly depends on pressure, leading to a pressure- 
dependent differential UCN production rate. A phenomenological expression for the multiphonon part of the scattering function 
s{X) describing UCN production has been derived from inelastic neutron scattering data. When combined with the production 
rate from single phonon processes this allows us to calculate the UCN production for any incident neutron flux. For calculations 
of the UCN production from single phonon processes we propose to use the values S* — 0.118(8) at saturated vapour pressure 
and 5** = 0.066(6) at 20 bar. As an example we will calculate the expected UCN production rate at the cold neutron beam for 
fundamental physics PFlb at the Institut Laue Langevin. We conclude that UCN production in superfluid helium under pressure 
is not attractive. 
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►v>( \ Introduction 

5— 1 ' Current projects to increase the density of ultra-cold ncu- 
5^ trons (UCN) available for physics experiments employ neu- 
tron converters of superfluid helium (He-II) I 'i .') and solid 
deuterium [4I5I6I7I8] . In these materials cold neutrons with 
energies in the order of meV may be scattered down to the 
neV energy range by the excitation of a single phonon or 
multiphonon in the converter. The inverse process is sup- 
pressed by the Boltzmann factor if the converter is kept at 
sufficiently low temperature. We consider UCN production 
in He-II, for which first publications have focused on the 
single phonon process, where incident 8.9 A neutrons are 
scattered down to the UCN energy range |9|10| . Later two 
publications [ll|12j have set out to calculate the expected 
UCN production rate from processes involving the emis- 
sion of two or more excitations but reached contradicting 
results. 

The purpose of this article is to provide on the one hand 
a physical-model-based extrapolation to short wavelengths 
of scattering data relevant for UCN production. On the 
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other hand we consider for the first time UCN production 
in He-II under pressure. This is motivated by the pressure 
dependence of the properties of He-II, in particular its dis- 
persion curve and density. Application of pressure to He- 
II increases the velocity of sound, such that the dispersion 
curves of He-II and of the free neutron cross at shorter neu- 
tron wavelength. For neutron beams from a neutron guide 
coupled to a liquid deuterium cold source, the differential 
flux density d$ /dA in the range 8 — 9 A normally increases 
for decreasing wavelength. Pressure also increases the den- 
sity of He-II. These two facts may lead to the expecta- 
tion that the single phonon UCN production rate increases 
with pressure. Furthermore, the multiphonon contributions 
might be favourably affected by application of pressure. 
This justifies to investigate whether pressure may provide 
an UCN source superior to a converter at saturated vapour 
pressure (SVP). 

The first part of this paper presents some general expres- 
sions for the UCN production rate in He-II as given in an 
internal note in 1982 by Pendlebury |13j and in Ref. [TT]. In 
the second part we consider UCN production due to mul- 
tiphonon scattering both for SVP and for 20 bar, using in- 
elastic scattering data for He-II measured at 0.5 K [14j. We 
give a formula approximating the contribution from multi- 



phonon processes to the UCN production rate. This might 
be useful to calculate, for any given spectrum of incident 
neutrons, the expected UCN production rate. As an exam- 
ple, this is done here for the cold neutron beam PFlb at 
the Institut Laue Langevin. 



UCN production rate, general expressions 
UCN production in superfluid helium is due to coher- 
ent inelastic scattering of the incident cold neutrons with 
energy E (and wavenumber k) down to energies E' ik'). 
The maximum final energy is defined by the wall potential 
of the converter vessel (252(2) neV for beryllium) with re- 
spect to the Fermi potential of He-II (18.5 neV at SVP): 
Vc = 233(2) neV. The production rate is given by 
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where A(j>/AE is the differential incident flux, N the helium 
number density, and Aa/AE' the differential cross section 
for inelastic neutron scattering. The latter is given by 
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where b is the neutron scattering length of ''He, huj — E — 
E' , q = k — k' , and S{q, hio) is the dynamic scattering 
function evaluated for values on the dispersion curve of the 
free neutron. Details of calculation can be found in Ref. 11 
where the general result is given in Eq. (9). We may write 
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where a = 1.34(2) barn [15] is the scattering cross section 
of *He, hkc ~ \/2mrt4, and 
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S{q,hLu)6{hLU~h^ky2mn)Auj (4) 



defines the UCN scattering function as function of the in- 
cident neutron wavelength A, making use oi q — k — 2-k /X, 
valid for k' k. We divide it into a single and a multi- 
phonon part, s(A) = si(A) + sn(A). 

The single phonon contribution can be approximated by 
si(A) = S*S{\* - A), where A* = ^ir/q* is the neutron 
wavelength at the intersection of the dispersion curves of 
the free neutron and the helium {q* = 0.706 for SVP), 
and S* = h /pg^jj S{q*,Huj)Auj denotes the intensity due to 
single phonon emission. Evaluation of the single-phonon 
UCN production rate yields 



where E* = h q /2mn, and the Jacobian factor 
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SVP 


bar 


A* [A] 


8.92(2) 


8.26(2) 


E* [meV] 


1.028 


1.20 


N [lO^^cm-S] 


2.1835 


2.5317 


13 


1.42(1) 


1.21(1) 


S* 


0.118(8) 


0.066(7) 


relative factor 


1 


0.41 



Table 1 

Factors relevant for single phonon production rate deduced from 
Ref. I16I17I181 . An increase in flux of a factor ~ 2.5, going from 8.9 A 
to 8.3 A is necessary to compensate for the loss in intensity. 
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accounts for the overlap of the two dispersion curves, En{q) 
and Enc (q) ■ For a beryllium-coated converter at SVP this 
results in 



PiiVc) = 4.97(38) • 10-8 A ^ 
cm dA 



(6) 



The differential multiphonon UCN production rate in the 
same units is given by 
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which depends explicitly on the form of sn(A) which is 
discussed below. 



UCN production rates from scattering data 
The scattering function Eq. ([4]) can be extracted from in- 
elastic scattering measurements. This has been done here 
for data previously measured by Andersen and colleagues, 
on the IN6 time-of-flight spectrometer at the Institut Laue 
Langevin, at SVP and at p = 20 bar. Parts of this data 
have already been published earlier in |14I18I19I20I21| . For 
the present purpose we have treated and rebinned all of the 
raw data to obtain a higher resolution in q and to cover the 
complete accessible g-range. Existing data points from pre- 
vious publications were used as reference for our analysis. 
A typical measurement of the dynamic scattering function 
S{q, huj) is shown in Figure[TJ 

The plane is divided into slices with width Aq = 0.2 A~^. 
Figure^] shows S{q, huj) for two representative values of q. 
From each such slice we extract one value of sii(A), with 
A fixed by the condition h/X = \/2m-ahL0. For this pur- 
pose the energy is binned with width Ahu) = 0.02 meV. 



The compiled data is presented in Table[2]of the appendix. 
For the short wavelength region below 4.5 A only one data 
point exists from Ref. [TS]. An accepted model for the ex- 
trapolation to short wavelengths is given by Family |22| . 
It was developed from low-order diagrams of the interpar- 




energy transfer ImeVl 



Fig. 1. Visualisation of S{q,hLj), 
p = 20 bar, measured by |14) . 
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Fig. 2. 5((?, ^iu;) for a) g = 0.90 A"! and b) ij = 0.95 A-^ The vertical 
lines indicate the energy of an incident neutron with E = cp' jlra^ 
that can be down-scattered to the UCN energy range. The width 
of the single phonon excitation is dominated by the finite resolution 
of the instrument. The roton-maxon (R+M) and two maxon (2M) 
resonances at higher energies are significantly lower in intensity. 

tide interactions in the high frequency Hmit and predicts 
S{q, hiS) cx g^w^''/^, which leads to sn(A) cx A'^ in our case. 

For a global characterisation of the multiphonon scat- 
tering function, we have found the following analytical ex- 
pression useful, 
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where i?MR is the sum of the maxon and roton energy. 

is the recoil energy of the He nucleus and 



2mHc 



/, p are fitting parameters which vary with pressure (see 
Table[3]for values). For A < 5 A it converges to the model 
of Family. In Figure[3] the correspondence of the model to 
the measured data is shown for SVP and 20 bar. 

Both the single and the multiphonon scattering functions 
change with pressure. The single phonon excitation moves 
to shorter wavelength, /3 and the value for S* decrease (see 
Eq. ([5|)). Therefore the UCN production by single phonon 
processes decreases as shown in Tabled! unless the cold 
neutron flux at 8.3 A is a factor of ~ 2.5 greater than at 
8.9 A to compensate for the loss in scattering intensity. 
The situation looks slightly different for the UCN pro- 
duction rate due to multiphonon scattering. With pres- 
sure the multiphonon excitations increase in intensity and 
the broad peak centre moves to shorter incident-neutron 
wavelengths (see Figure[3]). Integrating over Eq. ([7]) with 
the incident cold neutron flux of PFlb [23] at the Institut 
Laue Langevin, one can expect an increase in UCN produc- 
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Fig. 3. Multiphonon scattering function at SVP and 20 bar. The 
extrapolation to short wavelength of Korobkina et al. at SVP is 
linear in k, whereas the calculation of Schott et al. |12| is based on 
the static structure factor of superfluid helium. The data point (A) is 
taken from Ref. | 19l . The one-phonon peaks are indicated by vertical 
arrows: SVP (• ■ ■ ) and 20 bar ( — ). The data for 20 bar include the 
corrections proposed in Ref. |17| . 
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Fig. 4. Differential UCN production rate for SVP (— ) and 20 bar (- - -) 
at PFlb. The calculations have been performed for a critical wall 
potential of 252 neV. One can see clearly the shift of the single 
phonon peak to shorter wavelength, as well as the subtle increase 
from multiphonon processes. 

tion rate through multiphonon excitation. The differential 
production rate is shown in Figure[31 Integration over the 
UCN production curves for SVP and 20 bar gives in total 
(13.9±Q.9) cm-^s-i and (11.1 ±0.8) cm'^s'^ with a con- 
tribution of (9.5 ± 0.7) cm^^s-i and (5.8 ± 0.4) cm-^g-i 
from single phonon processes, respectively. 



Discussion 

We have shown how to calculate UCN production rates 
for different pressures based on inelastic neutron scattering 
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data. The single phonon UCN production rate is slightly 
higher than in Ref. pXi as we use for S* — 0.118(8), instead 
of S'* = 0.1. This higher value from Gibbs [18j (there called 
Z{Q)) is affirmed in a comparative analysis in Ref. [17]. 
A different approach by Yoshiki [24] , calculating the total 
cross section for UCN production from single phonon ex- 
citations in a beryllium coated converter, gives a slightly 
larger value. However, for his numerical calculation he did 
not take the Fermi potential of helium into account. 

A phenomenological expression is given to describe the 
global behaviour of the multiphonon scattering function. 
It converges for short incident wavelength to the model of 
Family. This model describes the short wavelength region 
better than the linear extrapolation in q used in Ref. [11 , 
which overestimates the multiphonon contribution. The 
calculations done by Schott and coworkers [12] predict a 
huge increase in the scattering function for short wave- 
lengths that appears to contradict the first moment sum 
rule [25j . Uncertainties remain in the region below 4 A 
as there is only one data point available 19] . Baker and 
coworkers jS^ have measured steadily decreasing production 
rates in the range of 6.5 A to 4.5 A. Therefore we assume to 
have only a small contribution to the scattering function for 
< 4 A, which we consider as being best represented by our 
expression. On the long wavelength end for A > A* there 
is no intensity, as there are no more elementary excitations 
in superfluid helium. 
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Conclusion 

UCN production in superfluid helium under pressure will 
not lead to a gain in production rate compared to SVP us- 
ing a cold neutron spectrum similar to the one of the cold 
beam PFlb. Only if the flux at 8.3 A exceeds by more than 
2.5 times that at 8.9 A, an increase in UCN production 
rate is expected, due to the increase of intensity of multi- 
phonon processes. However, one has to bear in mind that 
application of pressure requires a window for UCN extrac- 
tion. Whereas, an extraction window with associated severe 
UCN losses can be avoided for SVP [TT3] . 
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0.0580 ±0.0063 
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0.0359 ± 0.0062 





04984 ± 


00338 





91 


0.0410 ±0.0057 
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0.0626 ±0.0138 
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00394 
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0.0646 ± 0.0073 
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0.0746 ± 0.0088 





0574 ±0 


00384 
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05 


0.0628 ±0.0223 
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00425 
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07 


0.0851 ±0.0277 





06774 ± 


00419 


1 


10 


0.0674 ±0.0153 
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00442 
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0.0556 ±0.0058 
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0.0587 ±0.0061 





0556 ±0 


00431 
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0.0612 ±0.0062 





05191 ±0 


00555 
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0.0634 ± 0.0064 





04576 ± 
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0.0599 ±0.0065 





0427 ±0 
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0.0519 ±0.0058 





03515 ±0 
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0.0553 ± 0.0068 
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0.0488 ±0.0061 
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0.0500 ±0.0072 
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0.0473 ± 0.0064 
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0.0427 ±0.0063 
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0.0423 ± 0.0063 
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0.0486 ±0.0079 








1 


90 







014 ±0 
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Table 2 

Scattering function s(q) = h J S{q, huj)S(hw — h?q'^ /2mn)dLLj evalu- 
ated for UCN production at SVP and 20 bar. Blank spaces indicate 
that no data points are available for these wave vectors. 





SVP (T = 1.24 K) 


20 bar (T = 0.5 K) 


/ [lO-^^A-^meV-i] 
P 

Emr [meV-i] 


0.3 ±0.01 
4.46 ±0.26 
1.93 


0.6 ±0.02 
1.12±0.09 
1.94 



Table 3 

Emr from Ref. |18l and fitting parameters of Eq. ^ for SVP and 
20 bar, determined in a least— fit to the inelastic scattering data. 



